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We review our argument on the nature of the so-called meteorite constraint on the possible 
time-dependence of the fine-structure constant, emphasizing that dating meteorites at the 
present time is different in principle from searching directly for the traces in the past, as in the 
Oklo phenomenon and the QSO absorption lines. In the related literature, we still find some 
■ arguments not necessarily consistent with this difference to be taken properly into account. 

It does not immediately follow that any model-dependent approaches are useless in practice, 
though we cannot help suspecting that dating meteorites is no match for the Oklo and the 
' QSO in probing the time-variability of the fine-structure constant, at this moment. Some of 

. the relevance to the QSO data particularly in terms of the scalar field will be discussed. 

o , 

^ ■ 1 Introduction 

I Ever since Dirac's suggestion ^ that the gravitational "constant" G may not be truly constant, 

many measurements as well as theoretical analyses have been made on what had been accepted 
to be fundamental constants, as the experimental technique has steadily improved in the accuracy. 



p. 



D , Unfortunately his prediction G{t) ^ t ^ is already excluded by the recent observational findings, 

G/G = (0.2 ± 0.4) X IQ-^^y-^ and (-0.06 ± 0.2) x lO^^y^^ derived from the Viking Project and 
■ the analysis of the binary pulsars respectively. They are less than the predicted absolute value 

' I — ^0 ^1 0.7 X 10^^"y~^, by at least an order of magnitude, where to « 1.37 x 10^°y is the age 



of the Universe. An order-of-magnitude difference might look rather trivial compared with a huge 
number like tg. To appreciate how seriously it might affect the whole scenario, however, let us try 
to understand the above results by assuming an approximate equation G/G = —0.1 x i^"^, indicating 
G{t) t^"'^, which varies too slowly compared with the behavior ^ i^^, as will be shown. 

Dirac started with noticing a near coincidence between the two large numbers, today's value of 
the ratio of the electromagnetic to the gravitational coupling strengths of the electrons TZ{to) = 
a/(Gml) = 8TTa{Mp/mcf « 4.3 x lO''^ and to Ac « 3.3 x lO^^, where a = e'^/{4TThc) is the fine- 
structure constant, while Alp = {SttG/ (ch))^^/^ and nic are the Planck mass and the electron mass, 
respectively. Also tc — h/{mcC^) = 1.29 x lO^^^sec is one of the fundamental microscopic times. By 
assuming G{t) ^ t^^, he readily obtained a reasonably small value for TZ{tc) ~ 1.3 x 10**, sufficiently 
close to order unity, theoretically a more natural number than TZ{to). His "success" will be lost, 
however, if G{t) ~ t^'^-^, yielding TZ{tc) ~ 0.6 x 10^^, far from being anything close to one. 

These "failures" have, however, never deterred physicists from pursuing the idea of time-dependent 
constants. For example, TZ{to), one of the large numbers discussed above is no longer a serious issue 
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because we are now familiar with the renormahzation-group equations expressing the "hierarchy ratio" 
in terms of an exponential function, Mp/rric ~ exp(5/(2Q;)) with a comfortable choice b ^ 0.74. 

On the other hand, the current idea of "unification" has its own reason why we have to be concerned 
about possible time-dependent coupling constants among fields. In the fundamental theory in higher 
dimensions we may start with assuming truly constant coupling constants, which descend to the 
observed coupling constants in the real 4-dimensional world on being multiplied by a scalar field, like 
the one representing the size of internal space in the Kaluza-Klein theories or the "dilaton" in 

string theory [7] . Some of them might be identified as dark energy supposed to be responsible for 
the acceleration of the Universe |H|-E1- We expect that such a field evolves slowly with the cosmic 
time, and so do the observed coupling constants, in a way not necessarily like t^^. In this context, 
detecting time- variability of any of the coupling constants might even be viewed as an "evidence" of 
the presence of an underlying theory at the deeper level behind the world as we see it. 

Among various coupling constants and some of the mass ratios, the suspected time-variability 
of the fine-structure constant has been a special focus of the recent studies, because it is not only 
dimensionless but also allows an easier and closer access than others. Laboratory experiments belong 
obviously to this class of attempts, though they are designed to look at almost instantaneous rate of 
change at the present time [r5|. Measurements over much longer time-spans include the analysis of 
the Oklo phenomenon and the observation of the absorption lines from quasi-stellar objects (QSO). 

The Oklo phenomenon occurred at (1.95 ± 0.05) x 10^ years ago J^Ij providing an exceptionally 
stringent constraint ^7], either an upper bound Aa/a = (—0.8 ± 1.0) x 10^^ or a nonzero change 
(0.88 ±0.07) X 10~^ where Aa means the value of a at the Oklo time minus today's value. See the 
second of Ref. ^Sl for the slight revisions of the results in the first paper. The analysis of the neutron 
absorption process of ^''^Sm exploits an amplification effect due to (today's) very small resonance 
energy, Ero = 97.3meV, compared with the typical mass scale of nuclear phenomena of the order of 
MeV. Note that the timing corresponds to the fractional look-back time 5 = 1 — t/to « 0.142 with the 
age of the Universe to « 1.37 x lO^'^^y, which can also be viewed as the cosmological redshift z w 0.156, 
based on the now standard cosmological model with spatially flat 3-space with the cosmological 
constant parametrized by Ha — 0.7 and the Hubble parameter Hq = 72 km/sec/Mpc 'JW . 

Even with possible complications due to the strong interaction |^ |^ , it is quite unlikely that the 
upper bound turns out to be an order of magnitude larger than the Coulomb-only estimate, as was 
discussed in Appendix A of 9 , in which we only assumed the nearly complete cancellation between 
the energies from the Coulomb and the strong interactions for any value of Er unless a very special 
relation happens to be present between them. Hence more likely is the value below the Coulomb-only 
estimate. The conclusion is sufficiently general to survive including any specific mechanisms, like the 
time-dependent strange quark mass |22) . for example, in the strong-interaction contribution. 

It has been argued recently that the neutron flux might likely have been in the thermally non- 
equilibrium state, thus yielding a nonzero result Aa/a = (0.41 ''"'^'io.oe) ^ 10^^ at the Oklo time 
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Figure 1: The upper panel shows y — (Aa/a) x 10^ versus either the fractional look-back time 
s = 1 — t/to or the redshift z obtained from Ref. 25 (143 data-points), while the lower panel from 
Ref. PZI (case 1 in the first of the quoted references, 23 data-points). Notice an order-of-magnitude 
difference in the vertical scales shown on the left. The 2 data-points from Refs. |26[ 128) not shown 
here share nearly the same features as in the lower panel. Also shown are the Oklo and the so- 
called meteorite constraints. The error-bars on the former are invisibly small both in the vertical and 
horizontal scales, whereas those for the latter are controversial, being the main subject of this study. 

|23| . A more detailed analysis has been made basically in the same context 24 , reaching the upper 
bound (-0.05 ± 0.61) x lO"'^. As we find, the computed absorption cross section of samarium-149 
falls off toward the negatively large value of AEr = E,. — Ero for 2 billion years ago, not only for the 
Maxwellian distribution in Fig. 1 of |18j but also for the supposedly more realistic neutron spectrum 
in Fig. 3 of as well as in Fig. 12 of This makes it likely that the measured and the calculated 
cross sections meet together for AEr ^ — lOOmeV, hence Aa/a <, 0.9 x 10^^, as in our non-null result 
|18j . (No significant departure from the Maxwellian spectrum occurs for positive AEr.) Combining 
this with our argument above we find it unlikely that the Oklo constraint is large enough to be 
comparable with the current QSO result, as will be shown shortly, even if uncertainties are yet to 
be taken into account on the details of Reactor Zones (RZ) 10 and 13, from which we collected the 
samples to minimize the possible contamination from outside after the end of the reactor activity |18| . 

On the other hand, the QSO observations with the technique of the many- mult iplet method appear 
to indicate values at the level of |Aa/a| somewhere around 10~^ — 10^^, in the ranges 0.20 ^ s < 0.89, 
or 0.23 ^ z < 4.7 |2S]-|2HI- See Refs. [291 130) for reviewing the earlier attempts. The circumstance 
we are going to deal with may be illustrated in Fig. ^ in which we show the two sets of data. The 
143 data-points from Keck/HIRES ^ provides Aa/a = (-0.573 ± 0.113) x lO"'^ thus claiming a 
negatively nonzero change, while the other of 23 data-points from VLT-UVES |23 (case 1 in the first 
of the quoted references) reveals no such variation with (—0.06 ±0.06) x 10^^. Another group |26[I28| 
also from VLT-UVES favors a null result from their much more careful analysis of 2 data-points. It is 
even not clear if the contradictions come from some of the systematic uncertainties, or from the lack 
of sensitivity of the measurement implying that the values observed so far do not represent the true 
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effect supposed to be hidden at a deeper level. For the clarification of the last point of view, building 
larger telescopes is among future projects in this research area |28| . 

One of the common key issues is an apparent difference between the Oklo constraint and the 
reported QSO data by (far) more than an order of magnitude, corroborated particularly by our 
discussion outlined above. This may have something to do with the fact that the Oklo event lies 
outside the QSO range, as shown in Fig. ^ By exploiting this situation we may expect that Aa/a 
is not a monotonic function of time, like the constant values chosen by the simple weighted-mean 
estimate attempted in the conventional analyses [2S]-|2H1- It can even be oscillating in such a way 
that the Oklo time corresponds to one of zeros of the oscillation, also supported by an oscillatory 
behavior of the scalar field supposed to be responsible for the accelerating Universe 13-52 j a-s will 
be discussed briefly in the final Section. In this context, we could wish if we find another kind of 
phenomena which turn out to be useful to probe time-dependence of a. 

This is a reason why Dyson's analysis attempted decades ago still attracts attentions. He discussed 
a possible constraint obtained by examining the very long-lived beta decay ^^^Re — s-^^'^ Os in iron 
meteorites |3TI. By assuming a rather extreme circumstance, he obtained |Aq;/q;| < 10"^ 10~^ 
around the early time of the solar system. The bound is less stringent than those for the Oklo and 
the QSO which we are now interested in. This constraint, to be called the "meteorite constraint" 
in what follows, is nevertheless interesting because the crucial process might have occurred when the 
iron meteorites were formed about 4.6 billion years ago, also referred conveniently to the "meteorite 
time," which, as indicated in Fig. ^ lies between the Oklo time and most of the times for the QSO 
absorption lines. 

A possible importance of the rhenium decay appeared to be re-discovered when Olive et al. pointed 
out that the decay rate in question has been measured much more accurately than before, A ~ 
1.666 X 10~"'^^y~^, with the precision of 0.5% [32]. We believe that this approach needs more scrutiny 
by taking possible time-dependence of the decay rate into account correctly through an extension 
of the isochron analysis. We come eventually to a conceptually different approach, finding that the 
analysis tends to depending on details of how the decay rate varies with time, contrary to the model- 
independence claimed earlier. We expect corresponding reassessments to be made in the literature on 
the related subject, including Refs. [011^1001151] . 

We nevertheless try to apply the above extended analysis to find a way to reduce inherent uncer- 
tainties in any past studies. We present an explicit theoretical model, but obtaining the result not very 
much promising, in practice. These discussions appear to point to the view that dating meteorites 
may not provide a constraint as strong as the Oklo and the QSO, as was also suggested toward the 
end of III. A. 4 in [30] without much elaboration. This is probably a fair conclusion to be reached at 
this moment, unless an entirely new ingredient is added, as suggested briefly in the final Section. 

In Section 2, we begin with reviewing Dyson's classic work on the decay of rhenium-187. We then 
proceed to discuss the traditional isochron analysis in Section 3, based on which Olive et al. claimed a 
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new result nearly as stringent as the Oklo constraint. This inspired us to extend the isochron technique 
in Section 4, by allowing the decay rate X(t) to be time-dependent, due to the time-dependent fine- 
structure constant, leading us to the conclusion that the crucial role is played by the time-averaged 
decay rate A. The need to differentiate these two different kinds of decay rate is strongly emphasized. 
In Section 5, we first re-examine Dyson's analysis in the new light of the extended version of the 
isochron technique. We then propose a theoretical model which is expected to derive the change of 
A from that of A(t) based on an explicit but somewhat artificial mechanism, by fully exploiting the 
observational uncertainties still present in today's isochron analysis. The model turns out to select a 
time-range around the meteorite time, but with the resulting Aa/a short of being unique. In the final 
Section 6, we discuss the relevance to the Oklo constraint and the analyses of the QSO absorption 
lines, together with the related remarks on the scalar-tensor theory. 

2 Dyson's argument — I 

The exceptionally rare decay of rhenium-187 is a consequence of an exceptionally small mass difference 
of its beta-decay into osmium-187, Q — 2.5 keV, compared with the estimated Coulomb energy 
difference Qc — —15.8 MeV (according to the mass formula used by Dyson most of which 

happens to be canceled by the energy due to the strong interaction. In the Coulomb-only estimate, 
we make use of the relation Qc c>c a to derive 

a Qc Qc Q Q 
where AQ = Q{t) — Q{to), like the same convention for Aa. 
We find 

A oc QP, (2) 

with p w 2.84 according to ,31j (though p — 3 is used in [2] leaving most of the results practically 
unaffected). We then obtain 

^ + %^V-1, (3) 



A \ Q a 

where A is chosen to be today's value of the decay rate Ao ~ 1.67 x lO^^-'^y 
In the linear approximation with 



we use to find 

Aa _„ ,„_.AA 



- 5.6x10"''^^, (5) 
a A 



_1.78xl04^. (6) 



A 

The occurrence of a small (large) number like 10^^ (10^) on the right-hand sides represents an am- 
plification effect, as was emphasized in the previous Section on the Oklo constraint. 
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Figure 2: The decay rate is A2 in the time-interval between ti and ^2- 

As Dyson noticed in JSJ, if Aa/a is negative and sufficiently large, A A/ A on the left-hand side 
can be so large positively that the decay rate at the time nearly comparable with 4.6 billion years 
ago should have prompted rhenium to decay too rapidly to leave enough amount of rhenium as we 
observe today on the Earth. Dyson required this possibility to be ruled out. Note that, without any 
such change, we have (^o — ii)Ao ~ 0.077, hence e~(*''-*i)->^o ~ 0.926, implying that most of rhenium 
present at ii, the time of formation of meteorites, is still present now, where to — ii ~ 4.6 x lO^y. 

Since a is supposed to be different from today's value only in the past, A can be much larger also 
in the remote past, probably toward the initial time ti. As in Fig. |3 we may simplify the situation by 
assuming that the decay rate reached a significantly large value A2 only during an interval, ti < t < t2, 
which was the earliest time-interval toward ti: 

t2-ti 



1 

10' 



(7) 



to — ti 

for example, giving ^2 — ~ 4.6 x lO^y. Suppose A2 « 100 Aq, then we would find the product as 
large as {t2 ~ ti)\2 « (100/10) x 0.077 « 0.77, and hence e"*^'^^ = 0.46, which is dangerously small to 
be avoided. We choose the above A2 as representing a safe upper bound 

AA 



A 



< 100. 



Substituting this number into the left-hand side of we find 



Ac 



< 6.5 X 10" 



a 
a 



or - < 1.4 X lO^^^y"^ 



(8) 



(9) 



This is essentially Dyson's result which we re-interpreted in terms of the time-span ti < t < t2, 
for the sake of later discussion. In fact we will show that the presence of the new result from dating 
undermines the approach described above. Notice also that with the large value ((SJ we fail to justify 
the linear approximation, resulting in —Aa/a ~ 5.6 x lO"'^, which happens to be nearly an order-of- 
magnitude larger than the value given by 0. 

3 Isochron analysis — traditional 



Meteorites are considered to have been formed nearly at the time when the solar system was created, 
as fragments of small planets. A single core of each planet was fragmented into a number of different 
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meteorites. The original planets contained rhenium and osmium, which were created much earlier as 
supernovae had exploded. Some of the mother planets of meteorites were melted partially or totally 
and materials of metal and silicates were separated in the planets. The portions of liquid metal 
were crystallized to form solid metal cores, mainly composed of iron thus creating various different 
iron meteorites. Different temperatures and pressures, and hence different chemical fractionations 
provided initial contents of rhenium, for example, which were different from meteorite to meteorite. 
Iron meteorites ejected from the same mother planet have similar chemical properties, hence are 
categorized into the same group. Because the trapped rhenium-187 in iron meteorites decayed slowly 
into osmium-187, isotope ratios of rhenium and osmium changed gradually with time. 

Let ti denote the time when all the meteorites in the same group were formed. In each of them, 
we define the ratio A^rc — (^^^Re)/(^^®Os) in terms of the abundance of the stable isotope ^*^0s. We 
find 

iVReW =Ae-^(*-*^), (10) 

where A for the initial value differs from meteorite to meteorite, while A is for the decay rate. The 
decay product reappears as the abundance of ^*^0s also defined as the ratio to ^^*0s: 

Nos{t)^A{l-e-^^*-'^^)+B, (11) 

where B represents the common amount having been present already aX t = ti. From H10|l and (|11() 
the coefficient A is eliminated giving 

Nos{t)^S{t)Nno{t)+B, (12) 

where the time-dependent S{t) is defined by 

S{t) = e^(*-*i) ~ 1, (13) 

which starts with zero, increasing as t — ti. 

The advantage of using the expression H12|l will be appreciated if we measure the abundances of 
^^^Re and ^*^0s and plotted them in a 2-dimensional plot as in Fig. O often called an isochron plot. 
If they are found to align along a single straight line with a common slope S{t) and the intercept B, 
we may justify that the meteorites were born at the same time, hence belonging to the same group. 
The obtained value of the slope at the present time determines the product X{to —ti). If the value of 
A is known, then we find the age to — ti. This is the principle of dating. 

Smoliar et al. reported the latest development of dating iron meteorites |32) . They applied a very 
careful analysis of the data obtained from the four groups (IIA,IIIA,IVA, and IVB), determining the 
slope with the precision of 0.5%. On the other hand, the ages of the angrite meteorites assumed to 
be formed at about the same time as the iron meteorites have been known with better precision by 
using the same technique applied to the U-Pb, and Mn-Cr isochron pairs, partly due to much larger 
decay rates. Accepting the absolute age to — ti — (4.5578 ± 0.0004) x lO^y determined in this way 
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Figure 3: 2-dimcnsional plot of the Re-Os isoehron. 

|35j . we then find A for the Re-Os decay to be A « 1.666 x 10~^"'^y~^, with the precision of 0.5%, as 
mentioned before. This might be denoted by 

AA 



A 



< 5 X 10 



-3 



(14) 



Ohve et al. seemed to have expected that this new result implied a new constraint on the time- 
variability of a f2TJ, by assuming that the left-hand side of H14() is the same as, or at least the upper 
bound of the largest possible time-change of X{t) due to the time-dependent a{t). They substituted 
(ITU) into © to find 

< 2.5 X 10"^ (15) 



a 

which turns out to be nearly as stringent as the Oklo constraint, though the same authors later 
admitted (34) that their assumption cannot be justified in practice, conceding that (|15|l is "over- 
restrictive." 

They also assumed that (|15ll represents a time- variation that occurred around ti. This intrigued 
us because the time is somewhat older than the Oklo time, lying toward the near-end of the QSO 
range. The argument raised, however, the question on how they chose this specific time ti out of the 
entire history of the meteorites. This issue will be discussed in the next Section. 



4 Isoehron analysis — extended 

As we learned, most of the researchers working on the isoehron analysis have taken it for granted 
that the decay rate A is a pure constant, having had none of the slightest idea on its possible time- 
dependence. Notice, however, that the relation (|10|) is a consequence of the differential equation 



dt 



(16) 



where we have been immediately allowed to extend a constant A to a time-dependent A(t). This 
equation permits an obvious solution 



iVRo(i) = Aexp 



- / \{t')dt' 
Jti 



(17) 
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It is an easy exercise to substitute (|17f) into (|16f) finding the former solves the latter. 
With t = <o, we also put lfT7|) into the form of (^UJ) 

where A is defined by a time- average of \{t): 



A 



1 



\{t')dt' 



(18) 



(19) 



*0 - tl Jtl 

In this way we have come to find that the decay rate which has been determined to the precision of 
0.5% before \62\ is now re-interpreted as the time-varying decay rate averaged over the entire history 
of the meteorites [55] . 

' < 5 X 10~^. (20) 



A 



The notion of A was discussed briefly in the final published version of Ref. pi] . 




Figure 4: Two random examples, a straight line and a sinusoidal curve, of the function \(t) for the 
same average A. 

For a given averaged A, we may think of many different ways of time-dependence of \{t), as 
illustrated in Fig. 01 The average, no matter how precisely it may be measured, has nothing to 
do with how the function behaves. It masks any variation during the entire history between ti and 
Iq. In this sense, the dating technique made at the present time is of different nature from looking 
directly for the past events, like probing the ancient Oklo remnants or measuring the wavelength of 
the absorption lines coming from very distant objects as in the QSO observation. 

Apart from the question in principle, however, two quantities, averaged integral and in-between 
behaviors may have some relationship in practice. This ought to be model-dependent. We intend to 
look into some details expecting to understand better the physical implications. 

Olive et al. noticed 34 that today's observed value of A covers A with the uncertainty of about 
(1.5 — 3.0)% 12^]. They applied the simplest linear approximation, basically the same as the dashed 
straight line illustrated in Fig. ^ thus expecting |AA/A| <, (0.015 — 0.03) at the other end, t ^ ti. 
Aside from a fundamental question on how the linear approximation is supported by any physical 
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ground, we point out that the time-range during which X(t) deviates from A significantly is hkely 
quite wide as wc find in Fig. 01 in contrast either to a narrow time-range of the Oklo phenomenon 
or to nearly negligible horizontal error bars of each QSO data-point. Even with certain elaboration 
suggested in |34| . it seems by no means obvious how to confine the timing of the a variation into a 
narrow and well-defined range around the meteorite time, as also illustrated in Fig. 3 of |37) or in the 
Figure on page 61 of [311 for the larger audience.. 

They also attempted to rely on the negative Aa/a claimed by one of the QSO groups We 
should be cautious, however, in view of the apparently contradicting reports by other groups |26l 1271 
I28| . It is, as we observe, this type of conflict that is supposed to be solved by the meteorite constraint. 
Also related to this issue, one of us (Y.F.) has even proposed to fit the QSO data in terms of an 
oscillating a{t), which appears favored by the accelerating Universe as well as the Oklo constraint 
accepted to be more solid than the meteorite constraint This may further contribute to 

weaken the constraint on the time of the a changes on the basis of the meteorite dating. 





A 
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tJ tg tg t 

Figure 5: A2 is made the largest possible under the condition of the average A and positivity of A(t). 
e can be chosen infinitesimally small. 

5 Consequences of the extended analysis 

We present examples to show how different consequences the preceding Section entails. 
5.1 Dyson's argument — II 

We revisit Dyson's historical argument, which is expected to provide a lesson for other future studies. 
Fig. 121 will be replaced by Fig. [SJ in which X{t) for t > t2 has been set to zero, except for a small 
time-interval of the width e toward the present time to. This narrow portion is inserted technically to 
ensure that today's value is (close to) A, but e can be chosen infinitesimal mathematically. We chose 
A = between t2 and to — e in order for A2 to be as large as possible keeping the averaged value of 
the function X{t) to be A under the natural condition A(i) > 0. 
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Requiring the average to equal A results in the condition, in the limit e ^ 0: 

it2-h)X2 = {to-ti)X. (21) 

The right-hand side, which can also be put into {to — ti)/X^^ « 0.08 <C 1 , and so is the left-hand 
side, where we used ~ X^^ « 6 x 10^*^7. This implies that there was no time sufficiently long for 
rhenium to decay significantly before the time t2, no matter how large A2, hence how small a might 
have been. This illustrates how unlikely the naive argument, like the one due to Dyson, applies here, 
as far as the positivity condition of the decay rate is maintained. 

This critical comment would not have followed if A^^ were to happen unrealistically to be nearly 
comparable or shorter than to — ti, the age of the meteorites, because then (t2 — ti)X2 >, 1, and hence 
Re must have decayed out nearly entirely. We also point out that the "peak" between ti and t2 can 
be anywhere inside the range ti and to- This is one of the general features that knowing only the 
averaged decay rate tends to mask any details of the history. 

5.2 A simplified model relating AA/A to AA/A 

In this subsection we try to apply the extended isochron formulation to find a theoretical model sup- 
posed to be realistic, still appealing to the simplified step-function approach. We first point out that 
the straight line in the isochron plot, as in Fig. |31 is an approximation after all. In fact the isochron 
parameters, slope S and the intercept B, for each group have been determined in such a way that 
they minimize the chi-squared for the data-points on the individual iron meteorites (samples) in the 
group, each with the error bar. In Table 1, borrowed partly from Tables 1 and 2 of Ref. we show 
the number of samples and the estimated ages, relative to the absolute age of the angrites, as well as 
the 2(7 errors for each group. 

Table 1: Estimated ages and the number of samples of the groups |32| . 



Group 


IIA 


IIIA 


IVA 


IVB 


Number of samples 


10 


12 


10 


9 


Age (lO^y) 


4537 


4558 


4464 


4525 


Error (lO^y) 


8 


12 


26 


29 



These errors may be interpreted to be purely statistical in any conventional approach, but we dare 
to assume that individual iron meteorites in a group were indeed formed at different times distributed 
in the range indicated by the last row of Table due to some physical processes. This implies that 
a usual straight line in the isochron plot is only an approximate representation of 10 or so slightly 
different straight lines for individual meteorites. 

We also imagine that the different age yields different A thanks to a specific way of the time- 
dependence of A, as will be shown shortly. We would then be able to relate AA/A to AA/A. Let us 
consider an individual meteorite belonging to a group g. Also assume a rectangular-type function for 
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the time-dependent A(t): 



m = 




for h<t<t2, 
otherwise, 



(22) 



where 5 can be positive or negative but subject to the condition \5\ <^ Xq. We may have chosen 
different steps ^2 and ^3, respectively, at t2 and ^3, but decided not to do so for simphcity. Further 
comphcations might be avoided by choosing ^2 — ^3 ^ " ^i- It then follows 



AA 
A 



S_ 



(23) 



Now we assume that one of the meteorites is formed at tgi, which may differ from meteorite to 
meteorite but is confined to the range 



tn 



Wg ^ tgl <tg+ Wg, 



(24) 



where tg and Wg are taken from tlic second and the third rows, respectively, of Table Q It might 
be convenient to call this range the "initial 5-range." We also add another simplifying assumption 

A 



s 



g Si 

^--1 ^ 



ta tg 



to t 



Figure 6: The change of a taking place between t2 and before the initial g-range between tg 
and t„ + 



Suppose first the change of X{t) had taken place before tg 



as illustrated in Fig. El Obviously 



the "hump" (or dent if (5 < 0) does not contribute to the integral which appears in (|17(l with t = <o- 
Suppose, however, the right-end t2 of the hump enters the initial g-region between tg — Wg and tg + Wg. 
It may happen that t2 > tgi. We thus expect a change of the integral: 

(■to ^ 




\{t')dt' 



= \{t2~tg,)5\. 



Since the integral on the left-hand side defines (to — tgi)X, we find the change of A as given by 

it2-tgl)\S\ 



|AA| = 



tn-t 



9I 



(25) 



(26) 
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for a particular meteorite formed at tgi. This also results in a change of fluctuation of A of the group 
arising from different ages of the individual meteorites also due to the time-dependence of X(t). For 
an approximate estimate of the effect for the group g, we replace t2 — tgi by Wg as far as ^2 happens 
to land in the initial 5-range. Also by replacing tgi in the denominator by tg ~ ti, we obtain 



AA 



A 



AA 



A 



(27) 



to — ti A to — ti 

where we have used (I23II with the use of Ao ~ A. 

Suppose then the hump of 5 covers the entire initial g-range. In addition to the integral between 
tg — Wg to tg + Wg, wc must also include the contribution coming from the range from tg + Wg to 
t2. The former contribution remains essentially the same as 1)27(1 as far as the order of magnitude is 
concerned, while the latter is common to all the meteorites in the group g, hence fails to differentiate 
individual meteorites from each other. We obtain basically the same result ((27|l . In this sense the 
same contribution is obtained also when the left-end ts of the hump lands inside the initial g-range. 
Finally no contribution appears to A A/ A if the hump is entirely to the right of the initial g-range. 

According to Tabled we have Wg ^ 8.0 C x 10^ y where the coefficient C falls into the range roughly 
from ~ 1 for IIA to ~ 3.6 for IVB. We thus estimate 

«1.7CxlO"^ with 1<C<3.6. (28) 

to — ti 

We find it further necessary to estimate what portion of the "observed" value 0.5% comes from 
the S effect. We introduce another parameter < 1, in terms of which the left-hand side of H27|) is 
expressed by Sry x 10"'^ due to H2U|) . We also substitute from ^ expressing AA/A in terms of Aa/a. 
In this way we finally obtain 



Ac 



1.7 X 10 



c 



(29) 



which applies if a happens to be different from today's value in a range that covers the time of the 
meteorite formation, with a certain width ^ Wg. 

Now with a tentative choice ( — 1, we find lAa/aj ~ x 10~^. If r/ is close to one, the right- 
hand side of H29|l appears too large compared with the currently reported QSO data, certainly not a 
promising result. If choosing 77 « 1 were to turn out to yield somewhat smaller estimate than the QSO 
result, we would have been ready to accept our model as providing a stronger constraint than the 
QSO. Unfortunately that never happens to be the case. We are then forced to assume ij much smaller 
than unity, implying, rather naturally, that a major part of the uncertainty comes from measurement 
errors, for example. But then the result will be inevitably less useful. 

According to our previously repeated criticism, the analysis in terms of the time-averaged decay 
rate fails generally to locate precisely when a{t) changed. But a specific way of constructing the 
model has allowed us to focus selectively upon a narrow range around the meteorite time for a{t) 
which can have changed over a much wider range, though the resulting Aa/a is constrained in a way 
quite different from many of the related references. We will also comment in the next Section briefly 
on using H29|) as a test of overall consistency of the analysis. 
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6 Relevance to the Oklo and the QSO 



As we showed in Section 1, studying the meteorite constraint was motivated to give a solution to 
an apparent conflict between the Oklo and the QSO. From this point of view it seems pertinent to 
explain briefly what our overall scenario can be like, with special attention to H29|l . 

Notice first that, as was emphasized in Section 1, the "error-bar" for the Oklo constraint is at the 
level of lO^'^ — 10^^ in units oi y = (Aa/a) x 10^, as used in Fig. ^ If we extrapolate the horizontal 
straight line y — —0.573 for Keck/HIRES down to s = 0.142, this single point contributes 
more than hundreds, too large for the entire fit including the 143 QSO data-points to be acceptable. 
It was argued that the scalar field which is responsible for the time-dependent a at Oklo on the Earth 
is different from the one for the far distant objects (3^1. We find this, however, unlikely to be the 
case partly because we expect the space-dependent component of the scalar field to show the nature 
of finite force-range j40j in connection with non-Newtonian gravity. Rather we prefer a curved fitting 
function. 

A reasonable assumption might be provided by a damped oscillation [111 . 

y{s) = a (e^^^-i) cos(i; - vi) - e"'' cos(t;i)) , (30) 

where v/s — Uokio/soUo = 2nT^^ with the coefficient vi determined by 

vi = tan"^ ((e"^^""'" - cos(i;okio)) / sin(t;okio)) , (31) 

due to which y{s) vanishes naturally at s = (today) and s = Sokio ~ 0.142, as a good approximation 
to a very small value. We admit that this near identification of the Oklo with a zero of oscillation is 
provisional, still believing that the oscillation is supported by the accelerating Universe. 

Accepting H30I) as a phcnomenological function, at this moment, the best fit for the data set 
in [23] is obtained [TT| for a = 0.020,6 = 5.5, and T = 1.352 with x^d = 1-015, to be compared 
with Xrd ~ 1-023 for the original weighted-mean fit. Then at s = 0.33, the meteorite time, we find 
y = —0.072 which is about 3 times as large as the meteorite constraint \y\ — 0.025 corresponding to 
H15|) . In view of the non-uniqueness of this constraint as discussed before, no further effort appears to 
be worth attempting for s between the Oklo and the QSO times. We also find that trying to lessen 
the calculated \y\ for this range of s proves increasingly difficult unless it departs from the observed 
flat behavior excessively toward the high-s end ^1 EHl • 

The same technique applied to the data set of |27| reveals a different feature. The original weighted- 
mean analysis gave nearly a null result, whereas the data is best fit ,11, by the nonzero function (|30|) 
with a ^ -0.050, b = 3.1 and T = 0.134 yielding Xrd = 0-53, even better than the value 0.95 in [77| . 
The oscillatory behavior, being quite different from the nearly flat distribution in the data of 25 , is 
already apparent but can be easily washed out by a 1-parameter fit in terms of a horizontal straight 
line. Our fit gives y — 0.166 at s = 0.33, reaching close to a peak of the oscillation, however, with 
nearby zeros at s = 0.28 and s = 0.36. The meteorite constraint seems unlikely to play any crucial 
role even in its less restricted form. 
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Having estimated Aa/a from the QSO fits, we may substitute the resuhs back to the left-hand 
side of ^ to evaluate t], finding (0.5-1.8) x 10"^ a^^j (1.2-4.3) x 10"^ foj. ^ ^nd E?], respectively. 
No unique conclusion can be likely drawn from these numbers. We still hope that future more detailed 
analysis of the same kind might play a decisive role in selecting the true time-dependence of a obtained 
from the QSO absorption lines. 

We finally add how (|30|l is related to the cosmological acceleration. We start with pointing out that 
the scalar-tensor theory provides a successful implementation of the scenario of a decaying cosmological 
constant [5], according to which the observed A{t), time-dependent cosmological "constant" in the 
Einstein conformal frame, behaves like ~ t^^ (in the reduced Planckian unit system), allowing us to 
understand the smallness by 120 orders of magnitude without appealing to an extreme fine-tuning 
of parameters. Notice to ^ lO^*' in units of the Planckian time. In order further to understand the 
nearly constant behavior of A rather than falling off like t~'^, we need another mechanism to keep the 
scalar field to stay nearly standstill. Probably the easiest way is to assume that the scalar field is 
trapped by a potential which naturally causes a damped oscillation to be identified with what we find 
to fit the QSO data. This view is consistent, as we interpret jl2) . with the result of a careful analysis 
in Refs. j^EHI) which also appears to be consistent with the behavior encountered in |27) . 

Now from a more theoretical point of view keeping a strong tie with string theory than the 
phenomenology-oriented approach of quintessence |41| . the scalar-tensor theory may not allow us 
to ignore the emergence of a local field responsible for non-Newtonian gravity which features an in- 
termediate force-range of macroscopic size and violation of Weak Equivalence Principle (WEP). This 
is also expected to free us from the long-standing constraint from the solar-system experiment. For 
more details on these crucial issues, see Refs. |51 142| . 

It should be truly remarkable and fascinating if these different kinds of phenomena are related to 
possible time-variability of the fine-structure constant via a scalar field. We want to see if the issue 
of the rhenium decay contributes to verify these interrelations. 
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